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Abstract 

We proposed the procedure of measuring the unknown state of the 
three-level system - the qutrit, which was realized as the arbitrary po- 
larization state of the single-mode biphoton field. This procedure is ac- 
complished for the set of the pure states of qutrits; this set is defined by 
the properties of SU(2) transformations, that are done by the polarization 
transformers. 



1 Introduction 

In the physics of quantum information those systems that can be completely 
described in terms of three orthogonal states, were called qutrits (q-trit). In the 
case of a pure state the wave function of the three-level system can be written 
as: 

* = c 1 |l)+p 2 |2) + c 3 |3) ) (1) 

where |1), |2), |3) - are the orthogonal basis states. The complex coefficients Cj 
are the amplitudes of the basis states \j) and satisfy the following normalizing 
condition 

X>f = l- (2) 

Decomposition (1) is the generalization of the definition of qubit when the 
dimension of the quantum system d = 3. 

Several ways are known how to experimentally realize the multilevel quantum 
optical systems. In one of them [1] the interferometric procedure of state prepa- 
ration is used, when the attenuated laser pulses are sent into the multi-armed 
interferometer. The number of arms is equal to the system's dimensionality. 
Identification of basis states is done either through the pulses delay (temporal 
basis) or by the presence of constructive interference in the certain arm of the 
interferometer that is put into the registration system (energy basis) . The other 
example is the optical field that consists of the pairs of correlated photons that 
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belong to the different polarization modes. The preparation of such fields and 
their unitary transformations are reviewed in [2,3]. Multi-level systems attract 
a great interest in quantum cryptography, where by using those systems one 
can achieve secrecy growing when the eavesdropper uses so called symmetric 
individual attacks [4-6]. 

The problem of the adequate measurement of the parameters of the quantum 
state is one of the major ones in quantum information science. With the deci- 
sion of this problem one can expect the realization of the information's output 
devices, the protocols of error correction, quantum repeaters, and other quan- 
tum communication devices. From the fundamental point of view the question 
of minimal set of measurements that is needed for the complete description of 
the system's state is also very important. Let us notice that in some cases it 
is not necessary to perform a complete set of measurements to define system's 
purity [7]. 

Of course, for different types of quantum states one also uses different types 
of measurement procedures. As an example for a squeezed state of light, meth- 
ods of homodyne tomography are developing [8] ; they, in principle allow restor- 
ing the density matrix of n-photon Fock states [9] . For the polarization-squeezed 
[10] and scalar [11] light the fluctuations of Stokes parameters are analyzed and 
the quasi-probability function is restored [12]. In the case of two-photon fields 
one registers the set of the fourth order field moments in different spatial and 
polarization modes [13]. Let us notice that in context of each experimental 
procedure, the a-priori information about the properties of the examined state 
plays an important role. 

2 Biphotons as qutrits 

This work is devoted to the optical realization of the protocol that allows restor- 
ing the density matrix of an unknown three-level system quantum state. The 
polarization state of the two-photon field that belongs to the single spatial and 
frequency mode is the object of our investigation. In [2] it was shown that the 
pure state of such a field can be written as: 

* = ci |2, 0) + ca |1, 1> + C3 |0, 2) . (3) 

The two-photon Fock states in two orthogonal; polarization modes H and V 
serve as the basis states. So, for example, the second term in (3) corresponds 
to the existence of one photon in modes H and V with | C2 1 2 probability. The 
vacuum component |0, 0) is not considered in (3), because when the field is 
registered by the method of coincidence of photopulses, the contribution of this 
component to the measured correlation functions is equal to zero. The imaginary 
parts of the complex coefficients Cj — \cj \ exp{iipj} , j = 1,2,3 are the phases of 
the basis states. The total phase of the wave function is incidental, and that's 
why one of the phases can be excluded from the consideration, giving us the 
relative phases. For example, <pi 2 = (fi — y>2 and y> 32 = <p3 — f2- 
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For the description of the polarization properties of the single- mode biphoton 
field, in [14] the so-called polarization or fourth-order coherency matrix was 
introduced: 



A 


D 


E 


D* 


C 


F 


E* 


p* 


B 



The elements of this matrix represent the normally ordered fourth moments 
of the field that can be written as 

A = (at 2 a 2 ) , B = (tf 2 b 2 ) , C = (0+6+06) , 
D = (aWab) , E = (a+ 2 6 2 ) , F = (aW) . 

Here, a+ = a' H , 6+ = a/ v , a = an, b = ay - are the operators of photon 
creation and annihilation in polarization modes H and V, respectively. It can 
be noticed that the diagonal components of K4 are real. They characterize 
the intensity fluctuations in parallel (A and B) or orthogonal (C) polarization 
modes. Non-diagonal elements D, E, F in general case are complex. 

Since the state of biphoton field can be fully described by the fourth moments 
of the field, the elements of K4 matrix can be derived through the components 
of the density matrix of biphoton field. As an example for the pure state (3), 
by definition [15] p mn — c m c* n and 

P11 = |ci| 2 = A/2, P22 = |c 2 | 2 = C, p 3 3 = |c 3 | 2 = B/2, (5) 



pi 2 = c lC * 2 = D*/V2, P13 = c lC * 3 = E* /2, ^3 = c 2 c* 3 = F*/V2. (6) 
The condition of state's purity 

P 2 =P- (7) 

And normalizing condition 

Sp( P ) = l (8) 

impose the certain constraints between the elements of K4 matrix. Thus, from 
(8) it follows that 

A + B + 2C = 2, (9) 

And condition (7) gives 

E* = ABC/DF, \F\ 2 =BC, \D\ 2 = C(2 - B - 2C). (10) 

For the mixed state, the definition of density matrix also includes the com- 
plementary averaging with the classical distribution function P by the possible 
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states of the system, where P satisfies P» = 1, and the components of density 

»=i 

matrix are: 



Pmn — c mC n - (H) 

3 Measurement of the state of qutrits 

The question rises - how many (real) parameters should be measured to charac- 
terize completely the unknown state of biphoton field? From the definition and 
properties of the density matrix, it follows that in the case of the pure state, 
the number of real parameters, that define the state of the system that has a 
dimension d, is equal to 2d — 2, and in the case of a mixed state is equal to 
d 2 — 1 . Correspondingly for the qutrits in the first case one needs to know four 
real numbers, in the second case - eight. Considering that in experiment one 
measures the unnormalized state's amplitudes and conditions (8,9) are needed 
to be checked every time after measuring all three diagonal elements of K± ma- 
trix, we obtain that in pure state five moments are needed to be measured and 
in the mixed state - nine. 

Before we go further into the discussion of suggested protocol, we notice that 
the measurement procedure always leads to the destruction of our state that is 
caused by its interaction with the classical measuring device, in our case - with 
the detector. So when we speak about the input state, we always have in mind 
that it is introduced by a large enough set of copies, and part of them can be 
destroyed by the measurement. The results of the measurements will be applied 
to the rest part of the ensemble; this procedure lies at the heart of the ensemble 
method of quantum measurements. 

In quantum optics as the measuring apparatus for the fourth moments of the 
field usually serves the Brown-Twiss scheme, which consists of the beam splitter 
with the photo detectors in the output ports (Fig.l). Polarization transforma- 
tions in each spatial mode is done with the help of retardation plates (A/2 A/4) 
and polarization filters (polarizers) . Let's examine the normally ordered fourth 

moments of the field i?i 2 (xi, Q\, X2> #2) ~ (b'\b'\b' ib' 2^ , that are registered in 
the scheme that is shown at Fig.l, with the given wave plates orientation and 
the fixed polarization that is transmitted by polarizers. Our goal will be the 
search of the minimal set of such moments i.e. parameters that are measured 
in the experiment, from which one can compose all the elements of K4 ma- 
trix. In this case the input field will be transformed by the wave plates and 
polarizers in such way that the registered moments can be derived through the 
components of K4 matrix. Polarization transformations that are done by the 
wave plates are unitary and the polarizer plays the role of the polarization filter 
that sets the polarization state registered by the detectors. This idea is based 
on the classical schemes in which the Stokes parameters are measured [16]. It 
was also used in [13] for the measurement of the polarization properties of the 
biphoton light in two spatial modes (so-called two-qubit case). Let's notice that 
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the choice of quarter- and half- wave plates as the polarization transformers 
evidently is not unique, it is dictated by the convenience (these wave plates are 
most commonly used in the polarization experiments) and by the clearance of 
the transformations. 

The consecutive action of the beam splitter, two wave plates and the po- 
laraizer that transmits the vertical polarization on the signal (idler) photon is 
described by the following matrix transformations: 

^ y ^ = GvG\/ 2 G\/ 4: Gbs ( & ) ' ( 12 ) 

where a and b - are the annihilation operators of the input state in two orthog- 
onal polarization modes H and V at the input, a' and b' - are the annihilation 
operators at the output of the transformers; the state vector is written in Jones 
representation. 

°« - ( l /V ~ 2 l/VS ) < 13 > 

is the matrix that describes the action of the non-polarizing beam splitter, 

ov- ( s ; ) (i4, 

is the matrix of the polarizer that transmits the vertical field component, 

- ( l r . ) <*) 

are the matrices of the wave plates. Here coefficients t and r are equal to 

t = cos 8 + i sin 8 cos 2a, ,..„, 
r = i sin 8 sin 2a, 

where 8 - is the optical thickness, and a - is the angle between optical axis 
and the vertical direction (V). For the quarter and half- wave plates 8\/4 = 
7r/4, ax/4 = X, S\/2 = 7r /2 «a/2 = 0j an d these coefficients can be rewrit- 
ten as: 



t A/4 = 1/>/2(1 + *cos2x), 
r A/4 = i/V2sin2x, 



(17) 



t x/2 = icos2e 
r \/2 = tsm20. 

It is clearly seen from the definition of matrix how one can measure its 
diagonal components or moments A, B and C . In the first case the optical axes 
of all wave plates are set vertically - along the direction of the transmission of 
polarizer P. In the second case the polarization in two arms must be rotated 
by 90°, that is achieved by setting X i = O°,0i = 45°, \2 = 0°,6» 2 = 45°. In 



5 



the third case the polarization is rotated only in one arm, and according to 
scheme's symmetry it is not important in which. The set is: Xi — 0°, $i = 45°, 
X2 = 0°,9 2 = 0°. The more complex transformations are needed when non- 
diagonal components of K\ matrix are measured. As an example let's look how 
the following setting of elements acts on the input state does. Let xi = 0°, 9\ = 
45°, X2 = 45°, 6*2 = 22.5°. It is not hard to calculate that in this case 

R 12 = (b'\b'\b\b'2) = 1/8 [(a t2 a 2 ) - (a^ab) - (aW) + ( a t 6 t a& )] = 
1/8 [A + C - 2ReD] . 

(19) 

The measured moment in this case contains the contributions from three 
elements of the coherency matrix. Two of them are real diagonal components 
A, C . The third one is the real part of the (complex) non-diagonal element D. 
This example shows that since we are not directly measuring the phases of the 
states <pi, (f2 or ^3, but its cosine and sine, then the number of the measurements 
needed is increasing. In the real experiment, which description is shown below, 
in each arm of the Brown- Twiss scheme we used simpler configuration than the 
set of two rotating wave plates and the fixed polarizer. We considered the fact 
that when measuring the moments of the fourth order, the transformation that 
is done by the half wave plate and the fixed polarizer is equivalent to the action 
of one polarizer, which orientation is given by the angle (5. The rotation angles 
of half wave plate and the polarizer (3 are bounded by the equation: 

P = -29. (20) 

In the Table 1 we show the values of the orientation angles of the quarter 
wave plates (xi,2) an d polarizers (fti^) in two arms versus the value of the 
corresponding measured moment. This table essentially serves as the protocol 
of the reconstruction of the input state of the field that is presented by the 
biphoton-qutrit. It can be seen that in general case, the number of required 
measurements is equal to nine. First seven measurements realize the protocol 
for the pure input state. Two additional measurements are necessary for the 
definition of the cosine (sine) values of the corresponding phases. Eighth and 
ninth lines of the table show how to find the real and imaginary parts of the 
complex moment , which in the case of the pure state, according to (10), is 
derived through the rest of the moments. Let's notice that in our protocol for 
the definition of each non-diagonal elements of K4 matrix, only three moments 
must be known, what the minimal number of measurement needed is apparently. 

4 Experiment 

The experimental setup is shown at Fig. 2. It can be conventionally separated 
in two blocks - the block of the preparation of the input state and the block 
of measurements. First block includes the cw Ar+ laser, operating at 351 nm, 
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Tabic 1: The protocol of measurement of the moment set, which form the 
coherency matrix. 
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with the output power of 120 mW, which serves as the pump for the non- 
linear lithium iodate crystal, where the process of biphoton generation goes on. 
This block also includes the system of adjustment mirrors, quartz wave plate, 
which orientation could be smoothly varied, and the interference filter with the 
central wavelength of 702 nm, and 5 nm FWHM. Two-photon states of light were 
generated by the process of spontaneous parametric down-conversion (SPDC) 
inside the non-linear crystal. We used type-I collinear, frequency degenerate 
SPDC. The wavelength of biphoton radiation was X s = Xi = 2X p = 702 ± 9nm. 
The polarization of both photons was vertical. In this case, right after the 
crystal, the biphoton field was in state 

* = 4 n |0,2) + |uac). (21) 

The quartz wave plate was used to transform this state to the one that is 
described by equation (3) (Fig. 2). It is known that all the transformations that 
are done with the polarization of biphotons by the retardation plates can be 
described by the unitary (3x3) matrices G [2]: 



(22) 



where 




r 2 



V2tr 

G = G(S,a)=\ -V2tr* \t\ 2 - \r\ 2 V2t*r |, (23) 

-V2t*r* t* 2 

and coefficients t and r were introduced by (16). 

Matrices (23) give us the irreducible presentation of SU(2) group with 3x3 
dimensionality in the space of the state vectors (3). Let's notice that one cannot 
realize an arbitrary polarization state of a biphoton field, by using wave plates 
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only. In general case such transformations, together with the space of a state 
vectors (3) form the three-dimensional unitary presentation of SU(3) group. 

The thickness of the setting wave plate was h— 824±lmkm, therefore the 
parameter 5 = ? [n — n e )h was fixed. The second parameter a was changing 
during the experiment that allowed us to set the state of a biphoton field that 
was given onto the input of the measurement block. It is clear that all the states 
that were prepared in such a way did not drive our biphoton field out of the 
pure states class. 

The measurement block consists of the Brown- Twiss scheme that is shown 
at Fig.l. In our experiments we used Pockel cells instead of the quarter wave 
plates. The use of the Pockel cells seemed preferable to us, because it allowed 
controlling the polarization transformation distantly, by applying the certain 
voltages on them. The spectral control of the biphoton field was realized with 
the help of the spectrograph. Pulses coming from the detectors were driven onto 
the standard coincidence scheme, which measured the number of coincidences 
rate that is proportional to correlators (19). 

The measurement procedure was as follows. For the certain orientation of 
the setting wave plate, we performed a set of measurements that is described 
in Table 1. Then we rotated the setting wave plate by a, what corresponded to 
the change of the input state, and performed the same set of measurements. 

The dependence of modulo squared of the three state's amplitude and two 
phases from the rotation angle of a setting wave plate is plotted on Fig. 3-6. 
Each experimental dot on the plot corresponds to the certain input state that is 
given by a setting wave plate. Solid lines - result of theoretical calculations by 
formulas (22,23). Since all three measured moments contributed to the phase's 
calculation (look at Table 1); the errors of three measurements were added and 
the precision of a corresponding measurement was low. The main source of 
errors is the low quality of the Pockel cells and as a consequence - the inad- 
equacy of the polarization transformations that are done with these elements. 
Definitely by using the retardation wave plates that are working in zero-order 
interference regime is the only way to overcome this problem. We notion the 
good correspondence of the calculations and the experimental results when mea- 
suring the modules of state's amplitudes. All errors that appear here are due 
to the errors in setting the correct polarizer orientation angle. 

The presentation of the measured field states in terms of complex coefficients 
Cj as it is done on Fig. 3-6, is possible only for the pure states. For the mixed 
states, one needs to measure all six moments, that form the K4 (4) matrix. 
The real and imaginary parts of the moments D and F are shown at Fig. 7, 8. 
The calculated and measured components of density matrix of the state that 
correspond to the orientation of the setting wave plate by a — 25°, are shown 
at Fig. 9. Of course, for the pure states, which case was realized in our experi- 
ment, these moments can be derived through amplitudes Cj. Moment was not 
measured in this case, because as it is shown above, it could be derived through 
the other moments. 
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5 Conclusion 



In this work, we proposed the procedure of measuring the unknown state of the 
three-level system - the qutrit, which was realized as the arbitrary polariza- 
tion state of the single-mode biphoton field. This procedure is experimentally 
realized for the set of the pure states of qutrits; this set is defined by the prop- 
erties of SU(2) transformations, that are done by the polarization transformers 
(retardation plates). 

However, it seems to be interesting to realize the complete protocol of den- 
sity matrix restoration both for the pure and for the mixed states of qutrits. 
These experiments are now in progress and their results will be published soon. 
Separately, the question of maximum-likelihood estimation of the experimental 
results to the most probable quantity of p [13,17] will be reviewed. 

This work was done by the financial support of RFBR (grant # 02-02-16664) 
and INTAS (01-2122). 
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Figure 1: Brown- Twiss scheme in which one can measure the moments that 
form the coherency matrix K4. BS - beam splitter, RP - retardation plates of 
quarter and half wavelength, that are characterized by the parameters (5^/2) <^a/4 
and #a/2i Xa/4j P are the polarizers that transmit the vertical polarization. 
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Figure 2: Experimental setup. 
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Figure 3: Dependence of the parameter of the input state |ci| ,from the rotation 
angle of the setting retardation plate 
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Figure 5: Dependence of the parameter of the input state | C3 1 ,from the rotation 
angle of the setting retardation plate 
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Figure 6: Dependence of the parameter of the input state ip% — <p3, tf2 — <Pi, from 
the rotation angle of the setting retardation plate 
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Figure 7: Dependencies of the real and imaginary parts of the moment F from 
the rotation angle of the setting retardation plate. 
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Figure 8: Dependencies of the real and imaginary parts of the moment D from 
the rotation angle of the setting retardation plate. 
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Figure 9: Density matrix that is measure for the input state given by the angle 
of the setting wave plate a = 25°. Theoretically calculated values are placed in 
brackets. 
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